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We discuss the dynamics of interacting atomic bright solitons and dark bubbles in bulk immiscible 
Bose-Einstein condensates. Coherent matter-wave clusters can be constructed using dark-bright 
pairs with appropriate phases. In two dimensions we describe novel types of matter- wave molecules 
without a scalar counterpart that can be seen as bound states of vector objects. 
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Introduction.- One of the most remarkable achieve- 
ments in quantum physics in the last decade was the ex- 
perimental realization of coherent atomic matter waves 
with Bose-Einstein condensates (BECs) of ultra-cold al- 
kaline atoms [l|. Because of the inter-atomic interactions, 
BECs support different types of quantum nonlinear co- 
herent excitations: dark [2], bright Q and gap atomic 
solitons [4[ , vortices, vortex rings and related structures 

, shock waves Q , different types of vector solitons Q , 
and others Many other nonlinear phenomena have 
been theoretically studied in BECs [9[ . 

A challenging open problem both in BEC experiments 
and in nonlinear science at large is the construction of 
complex molecule-like coherent matter- wave structures. 
A possible strategy, that we explore in this letter, is to use 
atomic solitons and bubbles as bricks to construct stable 
matter wave aggregates which display phase-dependent 
properties. Related problems have been considered in 
nonlinear optics in relation with the propagation of op- 
tical beams in media with saturable nonlinearities [10j, 
it being very difficult to construct even metastable long- 
living soliton clusters. Nonlocal interactions allow for the 
creation of soliton clusters but local interactions such as 



those present in ordinary BECs pose great difficulties [11] 
In the field of matter- waves, trains with a finite num- 
ber of quasi-one dimensional bright solitons have been 
observed to be robust and stable [3[ due to the presence 
of the trap [12j. However, the idea does not work for 
higher dimensions due to the blow-up phenomenon 
With repulsive nonlinearities, dark solitons always repel 
each other and cannot form bound states [13[ . 

In this letter we will show how multicomponent 
homonuclear Bose-Einstein condensates in the immisci- 
ble regime allow for the construction of robust solitonic 
molecules. These matter- wave clusters display phase- 
dependent properties due to their coherent nature. 

Physical system.- We will consider two-component 
homonuclear BECs with atoms in two different hyperfine 
states |1) and |2) [14[. We will work in the immiscible 



regime and consider droplets of atoms in component |2) 
to be phase separated from a component |1) assumed to 
have a much larger number of particles. 

1-D atomic soliton molecules.- BECs tightly confined 
along two transverse directions are quasi-one dimensional 
and ruled in the mean field limit by the equations TB|, [l6| 
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for j = 1,2. Without loss of generality we will take the 
adimensional chemical potential /i = 1, and Ai = 0. Im- 
miscibility implies that g\ 2 > gn#22 • The normalization 
for u 2 is given by J R \u 2 \ 2 = 2(022/^0)^/2, «22 and N 2 be- 
ing the s-wave scattering length and number of atoms in 
1 2). Finally = ^h/muj± is the length-scale in which 
the adimensional spatial units are measured. Since we 
are interested on the bulk dynamics as in ring-shaped 
condensates no longitudinal traps will be considered. 

Let us first study the case of equal interaction coeffi- 
cients, gij = 1, i, j = 1, 2, that is very close to the realistic 
situation e.g. in Rb [HI or Na [ItJ- The phenomena to be 
described in this paper are not dependent on this specific 
choice of parameters and indeed , la ter we will consider 
the effect of tunning interactions [18|. In the former case, 
Eqs. (pQ) have explicit dark-bright soliton solutions 



u\ = i^s'ma + y7^ cosa ^anh (k [x — q(t)]) , (2a) 
u 2 = 



^i(f> ^iQt ^ikx tan a 



seek (k[x - q(t)]) .(2b) 



where n = ^/ ficos 2 a + (7Vb/4) 2 — Nb /4 is the inverse 
of the wave packet length, £1 = n 2 (l — tan 2 a) — Ao, and 
the soliton center position is q(t) = Xq + ttzsina [15| . 

The repulsive nature of the interaction between dark 
solitons prevents the formation of bound states of dark 
solitons. Concerning bright solitons, their mutual forces 
in the absence of external effects depend on the phase dif- 
ferences, A(j) = |0i — 02 1, going from repulsive for A0 = 
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FIG. 1: Interaction of two vector dark-bright solitons (1,2) for 
different amplitudes of the bright component Nb,i = Nb,2 (a) 
0.01 (b) 0.02 (c) 0.05. Other parameter values are given in 
the text. We plot pseudocolor plots of \u2(x,t)\ 2 ("bright" 
component) on the range x G [-20, 20], t G [0, 2000]. 
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FIG. 4: Interaction of three vector dark-bright solitons for 
N B = 1, X ,i = -8, X ,2 = 0, X ,3 = 8 with </>i = 0, 
<j)2 — 2.28 , 03 = 0. Other parameters as in Fig. [TJ Shown 
are pseudocolor plots of (a) \ui(x,t)\ and (b) \u2(x,t)\ on the 
range x G [-30, 30], t G [0, 1000]. 




FIG. 2: Interaction of two vector dark-bright solitons for 
Nb = 0.05 and Acj) = 3. Other parameters as in 
Fig. [1] Shown are pseudocolor plots of (a) \ui(x,t)\ 2 and 
(b) \u 2 (x,t)\ 2 on the range x G [-20,20],t G [0,2000]. 

to attractive for Acf) = tt. There is a critical intermedi- 
ate regime for Acj) in which unstable bound states can be 
constructed. 

Thus, a vector object including a dark soliton in one 
component and a bright soliton in another component 
could lead to a stable bound state with a second vec- 
tor soliton of the same type when the coherent inter- 
actions between the dark and the "droplet-like" bright 
component have opposite directions leading to what can 
be considered as a coherent molecule. 

We have studied numerically the interaction between 
two (hereafter denoted 1 and 2) initially static {ql\ = 
OL2 = 0) vector solitons given by Eqs. (j2j). Other param- 
eters are \i = 1, A2 = and the initial positions are 
Xq 5 i = 3, Xq,2 = — 3. In order to have attractive interac- 
tions between the bright components we choose Ac/) = tt. 
Our results are summarized in Fig. [TJ 

Small droplets of atoms of the second species are not 
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FIG. 3: Interaction of three vector dark-bright solitons for 
Nb = 0.05, Xo,i = —6, Xo,2 = 0, Xo^ = 6. Other parameters 
as in Fig. [I] Shown are pseudocolor plots of (a) t)\ 2 and 

(b) \u 2 (x,t)\ 2 on the range x G [-20,20],* G [0,2000]. 



able to stop the outgoing motion induced by the repulsive 
interaction between the black solitons in |1). However, 
for larger numbers of atoms in |2) we obtain a bound 
state of the two vector solitons, with a threshold num- 
ber of atoms of about Nb = 0.18. This is a remarkable 
result leading to a bound pair of dark-bright solitons in 
a system where all the interactions are repulsive due to 
the effective attraction provided by the phase between 
the "droplets" of atoms in component |2). 

The formation of these bound states is sensitive to the 
specific phase difference between the bright droplets in 
1 2) as it is seen in Fig. [2j where a change of the phase 
difference from tt to 3 suffices to destabilize the system. 

However, the addition of more components leads to a 
very interesting feature: although the atoms in |2) fea- 
ture a less trivial dynamics, they are still able to sustain 
the multisoliton bound state provided the beating pe- 
riod is faster than the typical evolution times of the dark 
component. An example is shown in Fig. [3j 

Curiously, multisoliton states are much more robust to 
perturbations in the phase differences (see Fig. [4]). We 
have also constructed atomic soliton molecules with more 
than three atoms, e.g. five soliton bound states, etc. 

Two-dimensional molecules.- When passing to higher 
dimensions the phenomenology changes essentially due 
to the fact that the more natural building blocks are 
vortices for the first component hosting a "droplet" of 
the second component. In the scalar case, moving two- 
dimensional bound states were found in Ref. [19] and 
correspond to vortex pairs of opposite circulation and 
rarefaction pulses. A similar phenomenology arises in 
the two-component case in the miscible regime [2o|. As 
the velocity increases, the distance between vortices of 
opposite circulation decreases to zero. The solutions at 
even higher velocity are localized density perturbations 
without zeros of the wavefunction. In two dimensions 
the sequence of solutions terminates with solutions ap- 
proaching zero energy and momentum as the velocity U 
approaches the speed of sound. 

Due to the structure of the vortex velocity field it is 
not possible to use the interactions of the bright part 
to construct bound states of equally charged vortices. 
Thus, unless a simple extension of the previous ideas is 
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not possible we will see in what follows that a novel phe- 
nomenology arises very different from the classical one 
described in Ref. |19j that can be linked to the presence 
of bound states of solitary wave structures. 

In what follows we will systematically construct soli- 
tary waves moving with speed U along the x— direction 
in two-dimensional two-component condensates in the 
phase-separation regime as solutions of the coupled GP 
system (pQ), where d/dt — >> — Ud/dx, and d 2 /dx 2 — >> V 2 

iU lt = l^i + ^-\H 2 -^\H 2 )^i, (3a) 

iU lt = ^ v2 ^ + (A-a|^| 2 -|^ 2 | 2 )^ 2 (3b) 

together with the boundary conditions — )• 1,^2 — >• 
0, as |x| — >• oo. In the phase separation regime a = 
912/911 = 912/922 > 1. Here A = /i 2 //ii where /ii 
and /i2 are the dimensional chemical potentials of ipi 
and ip2- We solve numerically the discretized version 
of Eqs. (j3j) by a Newton- Raphson algorithm combined 
with a secant algorithm to find A for a given constraint 
on N2 = f \ip2\ 2 dx dy . We obtain a family of solutions 
characterized by the velocity of propagation U, energy, 
E and impulse p = (pi + P2, 0), given by 
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The resulting families of solutions are given on Fig. [5] for 
various choices of N2 together with the Jones-Roberts 
(JR) dispersion relations [19[ for one-component conden- 
sates. For a given speed U, solitary solutions with higher 
a have lower energy and higher impulse. Firstly, in con- 
trast with the JR solutions, there is a stationary solitary 
wave corresponding to the ground state of the system 
with all the mass of the second component forming a ra- 
dially symmetric "bubble" in the centre of the depleted 
first component. This complex has a nonzero energy E gs . 
As the velocity increases from zero, the bubble becomes 
oblate in the direction of the motion with the velocity 
field of the first component is that of a dipole as Fig. [6] 
illustrates. There is a point on the dispersion curve where 
the velocity reaches its maximum - the inflection point. 
As energy and momentum continue to increase, the ve- 
locity decreases and the solutions become pairs of vor- 
tices of opposite circulation in the first component with 
the second component condensing in the vortex cores. In 
general the bubble-like solutions, for small E can be seen 
as a bound state of a JR rarefaction pulse and a mo- 
bile "filling" of the second component. Fig. [5jb) shows 
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FIG. 5: (Color online) (a) Energy (E — E gs )- impulse (jp — 
Pi +P2) dispersion curve of the solitary wave solutions of Eqs. 
for a = 1.2 and various choices of N2. The upper solid 
line corresponds to N2 = 0, i.e. the JR dispersion curve, 
shown here for comparison. To ease comparison we subtract 
the ground state energy E gs in the plot. Dashed black line 
corresponds to N2 = 8 and the dash-dot line corresponds 
to N2 = 20. (b) Velocity as a function of momentum for 
N2 = 20, the red curve corresponding to the JR case. The 
solid line corresponds to the JR case. 
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FIG. 6: (Color online) A stationary bound state solution for 
U = 0.2, N 2 = 80. (a) |V>i(x,?/)| 2 (b) \^ 2 {x,y)\ 2 . Shown 
are pseudocolor plots of the densities, arrows indicating the 
direction of the local current density j = — Re(z'0V'0*) and 
for (a) a few contour lines (black solid lines). The spatial 
region shown in the plot is (x, y) G [-30, 30] x [-20, 20] 



that there is a maximum velocity for the propagation of 
these solutions well separated from the sound speed ve- 
locity what can be physically interpreted as a signature 
of the mass of the second component, the heavier being 
the second component the smaller being this maximum 
velocity. To study the collisional properties of these lo- 
calized structures we have simulated the evolution of two 
bubbles set on a colliding course. Initially the bubbles are 
separated by a large distance, so that individually they 
are accurately represented by the solutions we found. We 
have observed several possible outcomes of such collisions 
summarized in Fig. [71 Almost identical slow colliding 
bubbles may form a bound stationary state even when 
they collide with an offset. Bubbles moving with large ve- 
locities may scatter at it/ 2 angle resembling the collision 
of two pairs of vortices of opposite circulation. Almost 
elastic collisions between these structures were observed 
when the velocities or masses of the bubbles were very 
different. We have observed that a bound state is more 
likely to be formed when bubbles have similar phases of 
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FIG. 7: Outcomes of coherent bubble-droplet pairs collisions. 
Density snapshots of the second component are shown be- 
fore (left column) and after (right column) the collision. In 
the first row, identical incident bubbles with U = 0.1, N 2 = 
40, p = 6.36, E — 2.99 form a bound state. In the second 
row, for U\ — U 2 = 0.2 incident bubbles with the same 
speed but different sizes emulate a elastic collision (left bub- 
ble: iV 2 = 40, E = 5.6, p = 14.9, right bubble: N 2 = 50, E = 
7.7, p = 20.1). In the third row, identical bubbles (parameters 
being U = 0.2, N 2 = 40, E = 5.6, p = 14.9) emulate vortex 
recombination. In the fourth row, identical bubbles (same as 
those in third row) collide with an offset of 8 dimensionless 
units and form a bound state. 

the second component and move slowly. In all such colli- 
sions, a small fraction of the total mass is lost and carried 
away by sound waves. The outgoing bubbles are solitary 
wave solutions as verified by energy-impulse calculations. 

Conclusions.- We have discussed the dynamics of inter- 
acting bubble-droplet pairs in immiscible Bose-Einstein 
condensates. Coherent atomic soliton molecules made 
up of dark-bright soliton pairs with appropriate phases 
can be constructed in ID scenarios. In two-dimensions, 
we have found novel types of robust bubble-like solitons 
without a scalar counterpart that can be seen as a bound 
state of two vector objects. Our ideas can be used to con- 
struct coherent atomic molecules made up of solitons. 
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